Complete organelle genomes of Korean fir, Abies koreana and phylogenomics of the gymnosperm genus Abies using nuclear and cytoplasmic DNA sequence data.

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus.

Molecular regulatory mechanisms of staminate strobilus development and dehiscence in Torreya grandis.

Gymnosperms are mostly dioecious, and their staminate strobili undergo a longer developmental period than those of angiosperms. However, the underlying molecular mechanisms remain unclear. This study aimed to identify key genes and pathways involved in staminate strobilus development and dehiscence in Torreya grandis. Through weighted gene co-expression network analysis (WGCNA), we identified fast elongation-related genes enriched in carbon metabolism and auxin signal transduction, whereas dehiscence-related genes were abundant in alpha-linolenic acid metabolism and the phenylpropanoid pathway. Based on WGCNA, we also identified PHYTOCHROME-INTERACTING FACTOR4 (TgPIF4) as a potential regulator for fast elongation of staminate strobilus and 2 WRKY proteins (TgWRKY3 and TgWRKY31) as potential regulators for staminate strobilus dehiscence. Multiple protein-DNA interaction analyses showed that TgPIF4 directly activates the expression of TRANSPORT INHIBITOR RESPONSE2 (TgTIR2) and NADP-MALIC ENZYME (TgNADP-ME). Overexpression of TgPIF4 significantly promoted staminate strobilus elongation by elevating auxin signal transduction and pyruvate content. TgWRKY3 and TgWRKY31 bind to the promoters of the lignin biosynthesis gene PHENYLALANINE AMMONIA-LYASE (TgPAL) and jasmonic acid metabolism gene JASMONATE O-METHYLTRANSFERASE (TgJMT), respectively, and directly activate their transcription. Overexpression of TgWRKY3 and TgWRKY31 in the staminate strobilus led to early dehiscence, accompanied by increased lignin and methyl jasmonate levels, respectively. Collectively, our findings offer a perspective for understanding the growth of staminate strobili in gymnosperms.

Large-Scale Identification and Characterization Analysis of VQ Family Genes in Plants, Especially Gymnosperms.

VQ motif-containing (VQ) proteins are a class of transcription regulatory cofactors widely present in plants, playing crucial roles in growth and development, stress response, and defense. Although there have been some reports on the member identification and functional research of VQ genes in some plants, there is still a lack of large-scale identification and clear graphical presentation of their basic characterization information to help us to better understand this family. Especially in gymnosperms, the VQ family genes and their evolutionary relationships have not yet been reported. In this study, we systematically identified 2469 VQ genes from 56 plant species, including bryophytes, gymnosperms, and angiosperms, and analyzed their molecular and evolutionary features. We found that amino acids are only highly conserved in the VQ domain, while other positions are relatively variable; most VQ genes encode relatively small proteins and do not have introns. The GC content in Poaceae plants is the highest (up to 70%); these VQ proteins can be divided into nine subgroups. In particular, we analyzed the molecular characteristics, chromosome distribution, duplication events, and expression levels of VQ genes in three gymnosperms: Ginkgo biloba, Taxus chinensis, and Pinus tabuliformis. In gymnosperms, VQ genes are classified into 11 groups, with highly similar motifs in each group; most VQ proteins have less than 300 amino acids and are predicted to be located in nucleus. Tandem duplication is an important driving force for the expansion of the VQ gene family, and the evolutionary processes of most VQ genes and duplication events are relatively independent; some candidate VQ genes are preliminarily screened, and they are likely to be involved in plant growth and stress and defense responses. These results provide detailed information and powerful references for further understanding and utilizing the VQ family genes in various plants.

Model Species to Investigate the Origin of Flowers.

The angiosperms, or flowering plants, arose at least 135 million years ago (Ma) and rapidly diversified to form over 300,000 species alive today. This group appears, however, to have separated from its closest living relatives, the extant gymnosperms, much earlier: over 300 Ma. Representatives of basally-diverging angiosperm lineages are of key importance to studies aimed at reconstructing the most recent common ancestor of living angiosperms, including its morphological, anatomical, eco-physiological and molecular aspects. Furthermore, evo-devo comparisons of angiosperms with living gymnosperms may help to determine how the many novel aspects of angiosperms, including those of the flower, first came about. This chapter reviews literature on the origin of angiosperms and focusses on basally-diverging angiosperms and gymnosperms that show advantages as potential experimental models, reviewing information and protocols for the use of these species in an evo-devo context. The final section suggests a means by which data from living and fossil groups could be integrated to better elucidate evolutionary events that took place on the long stem-lineage that apparently preceded the radiation of living angiosperms.

Context and Mutation in Gymnosperm Chloroplast DNA.

Mutations and subsequent repair processes are known to be strongly context-dependent in the flowering-plant chloroplast genome. At least six flanking bases, three on each side, can have an influence on the relative rates of different types of mutation at any given site. In this analysis, examine context and substitution at noncoding and fourfold degenerate coding sites in gymnosperm DNA. The sequences are analyzed in sets of three, allowing the inference of the substitution direction and the generation of context-dependent rate matrices. The size of the dataset limits the analysis to the tetranucleotide context of the sites, but the evidence shows that there are significant contextual effects, with patterns that are similar to those observed in angiosperms. These effects most likely represent an influence on the underlying mutation/repair dynamics. The data extend the plastome lineages that feature very complex patterns of mutation, which can have significant effects on the evolutionary dynamics of the chloroplast genome.

Exploring the Seasonal Dynamics and Molecular Mechanism of Wood Formation in Gymnosperm Trees.

Forests, comprising 31% of the Earth's surface, play pivotal roles in regulating the carbon, water, and energy cycles. Despite being far less diverse than angiosperms, gymnosperms account for over 50% of the global woody biomass production. To sustain growth and development, gymnosperms have evolved the capacity to sense and respond to cyclical environmental signals, such as changes in photoperiod and seasonal temperature, which initiate growth (spring and summer) and dormancy (fall and winter). Cambium, the lateral meristem responsible for wood formation, is reactivated through a complex interplay among hormonal, genetic, and epigenetic factors. Temperature signals perceived in early spring induce the synthesis of several phytohormones, including auxins, cytokinins, and gibberellins, which in turn reactivate cambium cells. Additionally, microRNA-mediated genetic and epigenetic pathways modulate cambial function. As a result, the cambium becomes active during the summer, resulting in active secondary xylem (i.e., wood) production, and starts to become inactive in autumn. This review summarizes and discusses recent findings regarding the climatic, hormonal, genetic, and epigenetic regulation of wood formation in gymnosperm trees (i.e., conifers) in response to seasonal changes.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis.

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C18 Δ9-elongase and a C20 Δ5-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ5-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

Evolutionary analysis of KED-rich proteins in plants.

During the course of evolution, organisms have developed genetic mechanisms in response to various environmental stresses including wounding from mechanical damage or herbivory-caused injury. A previous study of wounding response in the plant tobacco identified a unique wound-induced gene, aptly named KED due to its coding for a protein that has an unusually high content of amino acids lysine (K), glutamic acid (E) and aspartic acid (D). However, by far little is known about this intriguing gene. In this study, we investigated the evolutionary aspects of the KED-rich coding genes. We found that a consistent pattern of wound-induced KED gene expression is maintained across representative species of angiosperm and gymnosperm. KED genes can be identified in species from all groups of land plants (Embryophyta). All the KED proteins from vascular plants (Tracheophyta) including angiosperm, gymnosperm, fern and lycophyte share a conserved 19-amino acid domain near the C-terminus, whereas bryophytes (moss, liverwort and hornwort) possess KED-rich, multi-direct-repeat sequences that are distinct from the vascular plant KEDs. We detected KED-rich sequences in Charophyta species but not in Chlorophyta wherever genome sequences are available. Our studies suggest diverse and complex evolution pathways for land plant KED genes. Vascular plant KEDs exhibit high evolutionary conservation, implicating their shared function in response to wounding stress. The extraordinary enrichment of amino acids K, E and D in these groups of distinct and widely distributed proteins may reflect the structural and functional requirement for these three residues during some 600 million years of land plant evolution.

Complete mitochondrial genome of Thuja sutchuenensis and its implications on evolutionary analysis of complex mitogenome architecture in Cupressaceae.

BACKGROUND: The complex physical structure and abundant repeat sequences make it difficult to assemble the mitogenomes of seed plants, especially gymnosperms. Only approximately 33 mitogenomes of gymnosperms have been reported. However, as the most widely distributed and the second largest family among gymnosperms, Cupressaceae has only six assembled mitogenomes, including five draft mitogenomes and one complete mitogenome, which has greatly hindered the understanding of mitogenome evolution within this large family, even gymnosperms. RESULTS: In this study, we assembled and validated the complete mitogenome of Thuja sutchuenensis, with a size of 2.4 Mb. Multiple sequence units constituted its complex structure, which can be reduced to three linear contigs and one small circular contig. The analysis of repeat sequences indicated that the numbers of simple sequence repeats increased during the evolutionary history of gymnosperms, and the mitogenome of Thuja sutchuenensis harboured abundant extra-long repeats (more than 5 kb). Additionally, the longest repeat sequence identified in these seven gymnosperms also came from the mitogenome of Thuja sutchuenensis, with a length of up to 47 kb. The analysis of colinear blocks and gene clusters both revealed that the orders of mitochondrial genes within gymnosperms was not conserved. The comparative analysis showed that only four tRNAs were shared by seven gymnosperms, namely, trnD-GUC, trnE-UUC, trnI-CAU and trnY-GUA. Furthermore, four genes have undergone potential positive selection in most gymnosperm species, namely, atp8, ccmB, mttB and sdh4. CONCLUSION: We successfully assembled the second complete mitogenome within Cupressaceae and verified that it consisted of multiple sequence units. Our study also indicated that abundant long repeats may contribute to the generation of the complex conformation of the mitogenome of Thuja sutchuenensis. The investigation of Thuja sutchuenensis's mitogenome in our study provides new insight into further understanding the complex mitogenome architecture within gymnosperms.

Mega-sized pericentromeric blocks of simple telomeric repeats and their variants reveal patterns of chromosome evolution in ancient Cycadales genomes.

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms.

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Evolution of complex genome architecture in gymnosperms.

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5-15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants.

Evolution of mitochondrial RNA editing in extant gymnosperms.

To unveil the evolution of mitochondrial RNA editing in gymnosperms, we characterized mitochondrial genomes (mitogenomes), plastid genomes, RNA editing sites, and pentatricopeptide repeat (PPR) proteins from 10 key taxa representing four of the five extant gymnosperm clades. The assembled mitogenomes vary in gene content due to massive gene losses in Gnetum and Conifer II clades. Mitochondrial gene expression levels also vary according to protein function, with the most highly expressed genes involved in the respiratory complex. We identified 9132 mitochondrial C-to-U editing sites, as well as 2846 P-class and 8530 PLS-class PPR proteins. Regains of editing sites were demonstrated in Conifer II rps3 transcripts whose corresponding mitogenomic sequences lack introns due to retroprocessing. Our analyses reveal that non-synonymous editing is efficient and results in more codons encoding hydrophobic amino acids. In contrast, synonymous editing, although performed with variable efficiency, can increase the number of U-ending codons that are preferentially utilized in gymnosperm mitochondria. The inferred loss-to-gain ratio of mitochondrial editing sites in gymnosperms is 2.1:1, of which losses of non-synonymous editing are mainly due to genomic C-to-T substitutions. However, such substitutions only explain a small fraction of synonymous editing site losses, indicating distinct evolutionary mechanisms. We show that gymnosperms have experienced multiple lineage-specific duplications in PLS-class PPR proteins. These duplications likely contribute to accumulated RNA editing sites, as a mechanistic correlation between RNA editing and PLS-class PPR proteins is statistically supported.

The Cycas genome and the early evolution of seed plants.

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgo form a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycas genome contains four homologues of the fitD gene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensis contains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensis genome provides an important new resource of broad utility for biologists.

A novel cupulate seed plant, Xadzigacalix quatsinoensis gen. et sp. nov., provides new insight into the Mesozoic radiation of gymnosperms.

PREMISE: Anatomically preserved evidence for a novel clade of gymnosperms emphasizes diversity of seed plants immediately prior to the appearance of angiosperm fossils in the paleontological record. METHODS: Cupulate seeds from the Early Cretaceous Apple Bay locality (Vancouver Island) are described from serial cellulose acetate peels and three-dimensional reconstruction. Phylogenetic context is assessed through the comparative analysis of gymnosperm seed producing fructifications and maximum parsimony analysis of a revised morphological data set for seed plant phylogeny. RESULTS: Xadzigacalix quatsinoensis gen. et sp. nov. is characterized by an orthotropous ovule with an elongated micropyle and complex integument, enclosed within a radial cupule. The micropylar canal is elongated; and the nucellus extends into the micropyle to seal the post pollination ovule. Except at the apex of the micropyle, the seed is completely enclosed by a parenchymatous cupule with ca. 20 axially elongated secretory ducts. The cupulate seed is produced upon a triangular woody stele, consisting of a parenchymatous pith surrounded by radially aligned tracheids. The stele produces three short terete traces that terminate within the base of the cupule as transfusion tissue at the seed chalaza. CONCLUSIONS: Organography, vascularization, nature of the integument and nucellus, and configuration of the micropylar canal distinguish Xadzigacalix quatsinoensis from all other gymnosperm clades. Cladistic analyses suggest the new plant may have affinities with gnetophytes or angiosperms. These results are complemented with a critical re-evaluation of ovulate structures for Mesozoic gymnosperms, providing new insight into plant diversity immediately antecedent to the explosive diversification of flowering plants.

Towards the plastome evolution and phylogeny of Cycas L. (Cycadaceae): molecular-morphology discordance and gene tree space analysis.

BACKGROUND: Plastid genomes (plastomes) present great potential in resolving multiscale phylogenetic relationship but few studies have focused on the influence of genetic characteristics of plastid genes, such as genetic variation and phylogenetic discordance, in resolving the phylogeny within a lineage. Here we examine plastome characteristics of Cycas L., the most diverse genus among extant cycads, and investigate the deep phylogenetic relationships within Cycas by sampling 47 plastomes representing all major clades from six sections. RESULTS: All Cycas plastomes shared consistent gene content and structure with only one gene loss detected in Philippine species C. wadei. Three novel plastome regions (psbA-matK, trnN-ndhF, chlL-trnN) were identified as containing the highest nucleotide variability. Molecular evolutionary analysis showed most of the plastid protein-coding genes have been under purifying selection except ndhB. Phylogenomic analyses that alternatively included concatenated and coalescent methods, both identified four clades but with conflicting topologies at shallow nodes. Specifically, we found three species-rich Cycas sections, namely Stangerioides, Indosinenses and Cycas, were not or only weakly supported as monophyly based on plastomic phylogeny. Tree space analyses based on different tree-inference methods both revealed three gene clusters, of which the cluster with moderate genetic properties showed the best congruence with the favored phylogeny. CONCLUSIONS: Our exploration in plastomic data for Cycas supports the idea that plastid protein-coding genes may exhibit discordance in phylogenetic signals. The incongruence between molecular phylogeny and morphological classification reported here may largely be attributed to the uniparental attribute of plastid, which cannot offer sufficient information to resolve the phylogeny. Contrasting to a previous consensus that genes with longer sequences and a higher proportion of variances are superior for phylogeny reconstruction, our result implies that the most effective phylogenetic signals could come from loci that own moderate variation, GC content, sequence length, and underwent modest selection.

A revised view on the evolution of glutamine synthetase isoenzymes in plants.

Glutamine synthetase (GS) is a key enzyme responsible for the incorporation of inorganic nitrogen in the form of ammonium into the amino acid glutamine. In plants, two groups of functional GS enzymes are found: eubacterial GSIIb (GLN2) and eukaryotic GSIIe (GLN1/GS). Only GLN1/GS genes are found in vascular plants, which suggests that they are involved in the final adaptation of plants to terrestrial life. The present phylogenetic study reclassifies the different GS genes of seed plants into three clusters: GS1a, GS1b and GS2. The presence of genes encoding GS2 has been expanded to Cycadopsida gymnosperms, which suggests the origin of this gene in a common ancestor of Cycadopsida, Ginkgoopsida and angiosperms. GS1a genes have been identified in all gymnosperms, basal angiosperms and some Magnoliidae species. Previous studies in conifers and the gene expression profiles obtained in ginkgo and magnolia in the present work could explain the absence of GS1a in more recent angiosperm species (e.g. monocots and eudicots) as a result of the redundant roles of GS1a and GS2 in photosynthetic cells. Altogether, the results provide a better understanding of the evolution of plant GS isoenzymes and their physiological roles, which is valuable for improving crop nitrogen use efficiency and productivity. This new view of GS evolution in plants, including a new cytosolic GS group (GS1a), has important functional implications in the context of plant metabolism adaptation to global changes.

The Chinese pine genome and methylome unveil key features of conifer evolution.

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

De novo transcriptome analysis provides insights into the salt tolerance of Podocarpus macrophyllus under salinity stress.

BACKGROUND: Soil salinization is causing ecosystem degradation and crop yield reduction worldwide, and elucidation of the mechanism of salt-tolerant plants to improve crop yield is highly significant. Podocarpus macrophyllus is an ancient gymnosperm species with a unique environmental adaptation strategy that may be attributed to its lengthy evolutionary process. The present study investigated the physiological and molecular responses of P. macrophyllus plants to salt stress by analyzing its photosynthetic system and antioxidant enzyme activity. We also analyzed the differentially expressed genes (DEGs) in P. macrophyllus under salt stress using RNA sequencing and de novo transcriptome assembly. RESULTS: Salt treatment significantly affected the photosynthetic system in P. macrophyllus seedlings, which decreased chlorophyll content, altered chloroplast ultrastructure, and reduced photosynthesis. The activities of antioxidant enzymes increased significantly following salt stress treatment. Transcriptome analysis showed that salt stress induced a large number of genes involved in multiple metabolic and biological regulation processes. The transcription levels of genes that mediate phytohormone transport or signaling were altered. K+ and Ca2+ transporter-encoding genes and the MYB transcription factor were upregulated under salt stress. However, the genes involved in cell wall biosynthesis and secondary metabolism were downregulated. CONCLUSION: Our research identified some important pathways and putative genes involved in salt tolerance in P. macrophyllus and provided clues for elucidating the mechanism of salt tolerance and the utilization of the salt tolerance genes of P. macrophyllus for crop improvement.

Evolution and structural variations in chloroplast tRNAs in gymnosperms.

BACKGROUND: Chloroplast transfer RNAs (tRNAs) can participate in various vital processes. Gymnosperms have important ecological and economic value, and they are the dominant species in forest ecosystems in the Northern Hemisphere. However, the evolution and structural changes in chloroplast tRNAs in gymnosperms remain largely unclear. RESULTS: In this study, we determined the nucleotide evolution, phylogenetic relationships, and structural variations in 1779 chloroplast tRNAs in gymnosperms. The numbers and types of tRNA genes present in the chloroplast genomes of different gymnosperms did not differ greatly, where the average number of tRNAs was 33 and the frequencies of occurrence for various types of tRNAs were generally consistent. Nearly half of the anticodons were absent. Molecular sequence variation analysis identified the conserved secondary structures of tRNAs. About a quarter of the tRNA genes were found to contain precoded 3' CCA tails. A few tRNAs have undergone novel structural changes that are closely related to their minimum free energy, and these structural changes affect the stability of the tRNAs. Phylogenetic analysis showed that tRNAs have evolved from multiple common ancestors. The transition rate was higher than the transversion rate in gymnosperm chloroplast tRNAs. More loss events than duplication events have occurred in gymnosperm chloroplast tRNAs during their evolutionary process. CONCLUSIONS: These findings provide novel insights into the molecular evolution and biological characteristics of chloroplast tRNAs in gymnosperms.